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'■^ : Abstract. The recently discovered transiting very hot Jupiter, HAT-P-7b, a planet detected by 

I ' the telescopes of HATNet, turned out to be among the ones subjected to the highest irradiation 

from the parent star. As known, the combination of photometric and spectroscopic data for 

such an object yields the stellar, orbital and planetary parameters. In order to best characterize 

^ this particular planet, we carried out a complex analysis based on a complete and simultaneous 

r^ Monte-Carlo solution using all available data. We included the discovery hght curves, partial 

CIh" follow-up hght curves, the radial velocity data, and we used the stellar evolution models to infer 

J,^ ' the stellar properties. 

t^ , This self-consistent way of modeling provides the most precise estimate of the a posteriori 

"^ ' distributions of all of the system parameters of interest, and avoids making assumptions on 
j^ \ the values and uncertainties of any of the internally derived variables describing the system. 
I ' , This analysis demonstrates that even partial hght curve information can be valuable. This may 
become very important for future discoveries of planets with longer periods - and therefore 
s^ \ longer transit durations - where the chance of observing a full event is small. 
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1. Introduction 

^^ The transiting cxtrasolar planet candidate HTR154-011 has been identified by the 
OO 



telescopes of the HATNet (HAT-7, HAT-8 and HAT-6, HAT-9) - see |Bakos et al.(2002) 
for a detailed description about the instrumentation. The reduction of the frames fol- 
lowed the standard frame calibration procedures, followed by star detection, astrome- 
try and aperture photometry which was supported by various detrending algorithms. 
;J! ' We search the light curves for box-shaped transit signals using the BLS algorithm 
5^ ' (Kovacs, Zucker, & Mazeh(2002)l): a periodic dip was detected in the light curve of 
2MASS 19285935+4758102. The initial "rejection mode" spectroscopy with CfA Dig- 



ital Speedometer (Latham(1992) ) at FLWO eliminated the possibility of an eclipsing 
binary star system. The analysis of Keck HIRES spectra resulted the radial velocity 
data, stellar atmospheric parameters and the bisector information (which ruled out a 
"false positive" hierarchical triple). We also carried out a partial photometry follow-up 



in the Sloan z-band with KeplerCam on the FLWO 1.2m telescope. See Pal et al.(2008) 
for more details about the observational details and data reduction. 

The basics of the characterization of transiting extrasolar planets are the following. 
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The period of the orbital motion can be very accurately determined because of the long 
timespan of the discovery (HATNet) observations. The mass of the planet is proportional 
to the amplitude of the RV curve, K. The radius of the planet relative to the star deter- 
mines the depth of the transit (which is resulted both by the follow-up light curves and 
the HATNet data) . The inclination of the orbit can be derived from the shape of the tran- 
sit light curve, i.e. the impact parameter is proportional to the cosine of the inclination. 
To calculate the absolute mass and radius, we need the same data or the host star too, 
therefore we need at least one luminosity indicator for a good stellar evolution model. 



Because of the lack of parallax, this indicator can be a/Ri,, see Sozzetti et al.(2007) 
The bottleneck in this case is the partial follow-up, i.e. the shape parameters of the 
light curve - namely the ingress/egress length relative to the total duration - cannot be 
derived directly to obtain a reasonable value of a/R^,. 

2. System parameters 

To obtain the planetary, stellar and orbital parameters, we used all available data: (i) 
The HATNet light curves with approximately 16,600 data points (moderate precision 
with an rms of 7mmag). The HATNet data are affected by the period P, epoch Tccntcr 
(strongly), the planetary radius p (moderately), duration of the transit T^ui — '^{C/R*)^^ 
(moderately), and the impact parameter b (slightly), (ii) Radial velocity data, high pre- 
cision relative to the amplitude. The RV data are affected strongly by the period P, 
epoch Tcontcn semi-amplitude K, and eccentricity k = ecosw, h = esinuj. (iii) The 
partial follow-up light curve, acquired shortly after the radial velocity measurements. 
The follow-up photometry data are affected by the period P (slightly), epoch Tccntcr 
(strongly) , planetary radius p (strongly) , the duration of the transit Tdur (strongly) , and 
the impact parameter b (strongly). 

Because there was only a couple of periods time difference between the radial velocity 
observations and the partial follow-up photometry, the transit center time could have 
been very well extrapolated, thus the total transit duration can also be accurately ob- 
tained. There is a strong coupling between the adjustable parameters, moreover, the 
contribution of each parameter in each type of data are different (see above), therefore to 
obtain a plausible set of system parameters, a joint fit should be performed, incorporating 
all available data simultaneously. 

3. Monte-Carlo modeling 

The basic system parameters (see previous section) can be determined using a Monte- 
Carlo algorithm (the method of refitting to synthetic data sets, MCMC), resulting an 
a posteriori distribution of the adjustable parameters. This method also automatically 
results the a posteriori distribution of derived quantities, such as a/i?* or the inclination. 
The luminosity indicator is the density of the star p^, which is related to the observable 
quantities as 

3^ faY 3K faY 1 ^ _ Sq 

^*^GP^\rJ 2PGsini\Rj R.-^" RJ ^^''^ 

where both po f^nd Eq are observables, therefore these also have an a posteriori dis- 



tribution resulted by the Monte-Carlo fit. The stellar evolution model (Yi et al.(2001)) 
results the radius of the star as the function of the density, surface temperature T and 
metallicity z, i.e. 

R,^R4p,,T,z). (3.2) 



Obtaining planetary parameters using MC analysis 
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Spectroscopy 

Analyze HIRES data outside I2 re- 
gion. 



Isochrone fit 

Using Y^ isochrones with a/R-,, 
as luminosity indicator, determine 
lVh,R^,, age, log 3* -^ Is log 5* 
changed? If yes, fix new log(?*, go 
to SME analysis; if no go to final 
planetary parameters. 



SME analysis 

Determine Toff, logg*, [Fe/H], 

vsmi. Using T^s, logg*, [Fe/H] 

— > determine limb darkening 

parameters. 
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LC + RV modeling 

Joint fit to KeplerCam + HAT- 
Net LC and RV data -^ determine 
a/Ri,, Rp/Ri,, impact parameter b, 
transit center Tc and period P. 
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RV analysis 

Independent fit based on the RV 
measurements — > determine -ft'*, e 
and u!. In our case, the eccentricity 
found consistent with zero. 



Final stellar and planetary pa- 
rameters 

Based on all available spectroscopic 
and photometric data, obtain self- 
consistent values for stellar plane- 
tary, and orbital parameters. 



Figure 1. Flowchart of the iterative Monte-Carlo analysis of the stellar, orbital and planetary 
parameters. The arrows indicate the "data flow" between the subsequent steps. See text for 
further details. 

Equation p.ip and equation (|3.2p have two unknowns, i?* and p* since T and z are also 
observables, known from the spectroscopy (and also have a MC distribution from the 
SME analysis). Therefore, this set of equations can be solved, which yields the a poste- 
riori distribution of the stellar parameters. The combination of the stellar parameters 
Ri,, M* with the previously fitted quantities results the a posteriori distribution of the 
planetary parameters in a straightforward way. A sanity check can be done by comparing 
the resulted stellar surface gravity (obtained directly from M* and i?*) with the value 
provided by the SME. If they differ, the SME analysis has to be redone by fixing the 
surface gravity and the whole procedure can be done repeatedly until convergence, see 
flowchart on Fig. 1. We note that the complete MC way of analysis results the derived 
parameters independently from the form of calculations, e.g. the a posteriori distribu- 
tion of the semimajor axis is exactly the same if it is calculated as a = Ri,{a/ Ri,) or 
a = [G{Mi, + Afp)]i/3[F/(27r)]2/3, altough i?* or a/Ri, themselves have definitely larger 

1/3 
individual uncertainties than the period or M* . The final stellar, planetary and orbital 

parameters for the HAT-P-7(b) system are summarized in Table 1. 



4. Summary 

The outlined method provides a straightforward way for characterizing planetary sys- 
tems since (i) it omits any kind of intermediate values which otherwise would distort the 
error estimations, (ii) automatically results the uncertainties and correlations of the sys- 
tem parameters, including the stellar, planetary and orbital parameters too, (iii) it can 
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Table 1. Stellar, orbital and planetary parameters of the HAT-P-7(b) system. 



Stellar parameter 


Value 


Source 


Tcff (K) 


6350 ± 80 


SME 


[Fe/H] 


+0.26 ± 0.08 


SME 


log 5, (cgs) 


4.07 ± 0.06 


Y2+LC±SME 


vsini (km s^^) 


3.8 ±0.5 


SME 


M. (Mq) 


1.47 ± 0.06 


Y^±LC±SME 


R* (Re) 


IMtlil 


y2±lc±sme 


Mv (magnitude) 


3.00 ±0.22 


Y2±LC±SME 


Distance (pc) 


320±^g 


y2±lc±sme 


Orbital/planetary parameter 


Value 


P (days) 


2.204730 ± 0.000004 


Tcenter (BJD) 


2, 453, 790.2593 ±0.0010 


a/R, 




A QC+0.28 


p = Rp/R* 




0.0763 ± 0.0010 


b = acosi/Rt 




o-37j:°:i^ 


i (degrees) 




85°7ttl 


K (ms^^) 




213.5 ± 1.9 


e 




0.003 ±0.012 


Mp (Mj) 




1 TTB+O-"^^ 

J--' 'O-0.049 


Rp {R.i) 




1 QCO+O.IOS 


a (AU) 




0.0377 ± 0.0005 



utilize even the partial light curves which otherwise would not provide enough informa- 
tion for further steps - these partial follow-up light curves become even more important 
when systems with longer periods will be discovered (i.e. the chance to observe a com- 
plete transit on a single night will be definitely smaller). The planet HAT-P-7b is also 
very interesting since (i) it has a very high surface temperature, somewhere between 
2200 — 2700 K, depending on the heat redistribution, (ii) it happens to fall in the field 
of view of the upcoming Kepler Mission, yielding an opportunity for further continuous 
observations. 
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